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Abstract. This paper considers population transfer between eigenstates of a finite 
quantum ladder controlled by a classical electric field. Using an appropriate change of 
variables, we show that this setting can be set in the framework of adiabatic passage, 
which is known to facilitate ensemble control of quantum systems. Building on this 
insight, we present a mathematical proof of robustness for a control protocol - chirped 
pulse - practiced by experimentalists to drive an ensemble of quantum systems from 
the ground state to the most excited state. We then propose new adiabatic control 
protocols using a single chirped and amplitude shaped pulse, to robustly perform any 
permutation of cigenstatc populations, on an ensemble of systems with badly known 
coupling strengths. Such adiabatic control protocols arc illustrated by simulations 
achieving all 24 permutations for a 4-level ladder. 



1. Introduction 

Population transfer from eigenstate k to eigenstate I of a quantum system refers to 
finding a control input such that the projection of final system state on eigenstate I of the 
free Hamiltonian has the same norm as the projection of initial system state on eigenstate 
k. Applications of population transfer range from population inversion [I], where k and 
I are lowest and highest energy eigenstates, to quantum information processing [2j [3j 0] , 
where logic gates would (selectively) permute the populations of several eigenstates. 
In many applications, including those mentioned, relative insensitivity to variations in 
system parameters is important for robustness issues. 

In the present paper, we show how control inputs designed on the basis of 
adiabatic passage can implement any given permutation of eigenstate populations for 
a finite anharmonic quantum ladder. The controls we use are chirped pulses [5] with 
appropriately modulated amplitudes and exploit the idea of eigenvalue crossing [6]. 
The ladder consists of a free Hamiltonian with approximately equidistant eigenvalues 
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and where the control input couples eigenstates associated to consecutive eigenvalues. 
A striking robustness feature is that our control fields must only satisfy a set of key 
properties and achieve population transfer independently of the values of dipole moments 
coupling consecutive levels of the ladder. This is a major difference with respect to early 
non-adiabatic approaches to molecular ladder dissociation using chirped pulses [5]. 

In this sense, we achieve a specific form of ensemble control. Ensemble control in 
its most general form wants a same input to drive an ensemble of systems, with different 
values of some parameter p, from given p-dependent initial state to given p-dependent 
final state [7J definition 1]. Currently, solutions to this general problem are essentially 
restricted to two-level systems, achieving approximate ensemble control in finite time 
and exact ensemble control in infinite time [3 [9] . They rely on accurate knowledge 
of laser-system coupling strengths and accurately tailored inputs, involving e.g. exact 
instantaneous "7r-amplitude-impulses" . In our setting, system parameters need not be 
exactly known and the input must only satisfy a few key properties. In turn, regarding 
initial-to-final-state transformations, we are limited to population permutations (with 
arbitrary relative phases between components of different eigenstates) that are constant 
as a function of system parameters. Driving an ensemble of 2-level systems from a 
common initial to a common final state has also been much studied in the NMR context, 
e.g. with geometric methods [8]. 

Adiabatic passage is a control strategy that builds on the adiabatic evolution 
property: A system state initially close to an eigenstate of a time-varying Hamiltonian 
H(t) approximately follows the time- varying eigenstate of H(t) if it varies slowly enough; 
the slower H(t) varies, the better the adiabatic approximation. A thorough formal 
study of adiabatic evolution can be found in [TQl CHI EE], on which we build the 
proofs of our results. Adiabatic evolution has been standard since the early days of 
quantum mechanics [13], e.g. when interpreting system evolution in terms of "avoided 
eigenvalue crossings". In a ladder control context, population inversion in two-level 
systems by a "chirped" pulse — where frequencies of a Gaussian laser pulse are spread 
out in time — is known by experimentalists and theoretically explained in the adiabatic 
framework [H]. This is the most basic case of our control, section [3] with N = 2. Many 
experimentalists have then focused on multiple-laser techniques, individually addressing 
pairwise couplings in an A r -level system; this includes stimulated Raman adiabatic 
passage (STIRAP), see e.g. [151 ESJ 02] • For iV-level ladder systems specifically, the 
possibility of population transfer from the lowest to the highest energy eigenstates with a 
single chirped laser pulse has been recognized and exploited in "adiabatic rapid passage" 
experiments [THJ [TJ HU [17] . An analysis of N- level adiabatic molecular dissociation with 
chirped pulses is given in [20] based on the Floquet representation. In the present 
paper we provide a simple mathematical proof of population inversion with avoided 
crossings (gap condition) based on Favard's Theorem [21] and on the roots of orthogonal 
polynomials [22], and extend the framework by adding amplitude control to perform not 
only population inversion but all different permutations of free Hamiltonian eigenstates. 

The paper is organized as follows. Section [3] gives the formal statement and section 
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the proof for iV-level population inversion with "adiabatic rapid passage" , actually 
proving how initial population of level k is finally transferred to level iV — k — 1 in 
adiabatic approximation. The key point for using adiabatic passage is a change of 
frame that depends on time- varying control input phase; it is detailed in section [2] after 
formal description of the ladder system. The proof then applies the standard "adiabatic 
theorem with spectral gap condition", where time- varying eigenvalues are shown to 
remain separated for all times. The inversion is insensitive to exact energy values of the 
individual levels in the ladder. Section H] proposes adiabatic control inputs to transfer 
population between two arbitrary eigenstates. It requires the control field to vanish at 
specific times which depend on (some) energy levels of the anharmonic ladder, such 
that we select a pair of time-varying eigenvalues to cross. System evolution is then 
ruled by the "adiabatic theorem without spectral gap condition". A complementary 
study of system behavior in the neighborhood of two crossing eigenvalues and valid for 
more general systems than ladder ones, can be found in [23] . We again provide a formal 
proof of the control's effect and highlight its ensemble/robustness features in section El 
Section [5] finally shows how any permutation of eigenstate populations can be achieved 
in this adiabatic passage framework. Each control protocol is illustrated by a simulation 
at the end of the corresponding section. 

Notation: We use the Dirac bra-ket notations: G denotes a complex vector, 
{ip\ = ' is its Hermitian transpose, and (.|.) : 

\if} 2 ) is the Hermitian scalar product. For z G C we note $l(z) its real part and z* 
its conjugate. Hn is the set of N x N Hermitian matrices, where N G N. We note / 
the N x N identity matrix. For any matrix A G C NxN , we denote its Frobenius (or 
Hilbert-Schmidt) norm ||A|| = Vtr A^ A where tr • denotes trace. For H G "Hat, it holds 

11-^11 = A? where A , . . . , Ajv-i are the (real) eigenvalues of H. For H G Hn an d 

A an eigenvalue of H, we denote P\ eHn the orthogonal projector on the eigenspace of 
H associated to the eigenvalue A. If H has M distinct eigenvalues {Ao, .., Aju-i}, with 
M < N, then H = ^2^=0 ^kP\ k is the spectral decomposition of H. If M = N, then 
H is called non degenerate and each P\ k is a rank-one projector. When M < N we say 
that H is degenerate; then some P\ k have rank larger than 1. 

S 1 denotes the unit circle equivalent to K modulo 2n. For J an interval of R, the 
derivative of a differentiable function / : J — > S 1 is a function from J to R. For all 
n G N, we denote C n (J,K) the set of n times continuously differentiable functions 
from J to K, where J is an interval of R and K is an interval of R or S 1 . A multi- 
component function is n times continuously differentiable, e.g. H(s) G C n (J, "H^), if all 
its components belong to C n (J,K). For / G C\J,K C R n ), we note f'(y) G C°(J,R n ) 
the value at y G J of the derivative of /. R>o is the set of strictly positive real numbers; 
we use analog notation with >, < or <. is the set of integers from a to b, both 
boundaries included. When writing Co, . . . , cjv_i G S we mean that Ck belongs to the set 
S for each k G N^" 1 . Infimum and supremum of a set are noted sup and inf respectively. 
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2. Problem setting 

2.1. Standard formulation 

Consider a quantum system with wavefunction G C N , = 1, N G N, whose 

dynamics is governed by the Schrodinger equation (with h = 1) 

ij t m)) = {H Q +u{t)H l )\m)- (i) 

The Hamiltonians H £ H^ and Hi £ H N respectively characterize free and control- 
induced evolution, u(t) being a real scalar control. In the present paper, we consider a 
quantum ladder for which the Hamiltonians, in the eigenbasis {|0), . . . , \N — 1)} of H , 
take the form 

N-l 

H = J2 k (u + A k )\k)(k\ (2) 

fc=0 
N-2 

Hi= ^2fjL k (\k)(k + l\ + \k + l)(k\) , (3) 

fc=0 

with ujq G R >0 ; A , . . . , Ajv_! G R; and fi , . . . ,fiN-2 £ R>o- We assume that system 
([I]) features two very different orders of magnitude, 

w |A fe | <C uj for all k and all t . (4) 

Physically, H is the free Hamiltonian of a quantum ladder with mean resonant 
frequency uo and anharmonicities We call eigenstates |0), . . . , \N — 1} of H Q the 
levels of the ladder. Hi is the dipole moment matrix and models couplings between 
consecutive eigenstates; it is therefore tridiagonal with zero diagonal elements, and can 
be taken real positive and symmetric without loss of generality. Condition (j4j) expresses 
that control amplitude is relatively weak and that the ladder is close to a harmonic 
one, i.e. eigenvalues of H associated to consecutive eigenstates are close to equidistant. 
This allows to exploit resonant transitions between all consecutive eigenstates with a 
control of carrier frequency ujq. We consider a typical such control with a small positive 
parameter e, 

u(t) = 2 3? (e ZU)at E(t)) , E(t) = A(et)e^ et) (5) 

with \\j t E(t)\\ <w , (6) 

where A(t) G R and 9{t) G S 1 for all t G R>o- Parameter e governs the rate of variations 
in the envelope A(et) and frequency -^^d(et) = 9'(ei) of E(t); we show in the next 
sections how taking e small allows to apply adiabatic passage properties. The slow but 
nonzero frequency variation is a key element for our control strategy. Physically, control 
fields like ([5]) are obtained e.g. by "shaping" a single laser pulse [24~| . 

The rotating wave approximation (RWA), standard in quantum systems modeling, 
consists in writing (TI]) with the change of variable \4>(t)) = (j2k=o e tkuJot \k) (k\ J \ip(t)) 
and neglecting fast oscillating terms, to keep only those that vary at frequencies oj . 
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It can be justified by averaging theory [25] thanks to inequalities flU),®- Within this 
approximation, \<f>) follows the dynamics 

ij t m)) = {H I +H I {t))\m (?) 



where 



N-X 



Hj =J2 kA k\k)(k\ 

k=0 
N-2 

Hj(t) = Vk(E(t)\k)(k + 1\ + E*(t)\k + l)(k\) . 



k=0 

2.2. Change of frame 

Hamiltonian Hj(t) contains a control field whose phase \9{si) varies on timescales 
of order one. The key idea to apply adiabatic passage to the iV-level system is an 
appropriate further change of frame on (0), such that all explicit time-dependence 
in the resulting dynamics involves timescales of order e. To this end, we extend the 
change of frame given in [TJ] Section 4.6] for the two-level case and define = 
^n-x e ki di£t ) |^^_ Dynamics ([J]) becomes 

ij t \m) = (H R Hst)) + A{et)H 1 ) |£(f)> (8) 
with u = 9', Hi given by (JHJ) and 

N-X 

H R (v) = k ( A k - v ) \ k )( k \ for all v e R . (9) 

k=0 

Define the propagator U e to be a time- dependent N by N unitary matrix such that 
the solution of © is given by |f(t)) = U e {t)\£{0)) for all t and for all |{(0)). Then U e 
follows the dynamics 

ie^-U £ {s) = H(s) U £ (s) , U £ (0) = I (10) 
as 

with H(s) = H R {u(s)) + A(s)H 1 (11) 

in the time scale s = et. In the following, we study system ([TO]) for s in the interval 
[0, 1] and with A(s) and w{s) as controls. Our goal is to achieve: 

(a) Adiabatic approximate eigenstate permutations: 

lim max \\ U £ (l)\k) (k\U £ -\a(k)) (a {k)\\\=0 (12) 

e->0+ keG 

for given G C N^" 1 and given permutation a of (0, . . . , iV — 1). 

(b) Ensemble control: a single control (A, u) achieves such eigenstate permutation on an 
ensemble of systems with different parameter values; the parameters are the dipole 
moments (/x , • • • , f^N-2) an d, in some cases, the anharmonicities (A , . . . , Ajv-i). 

(c) Robust control inputs: the above holds for any (A, u) that satisfy a set of key 
properties. 
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Remark 1. Writing ( TTDj) in terms of\k)(k\, the projector on eigenspace {/3 \k) : /3 G C} 7 
expresses that the goal is really population transfer, i.e. we allow U e (l)\k) « e* Xfe |o"(A;)) 
w«£/i arbitrary phases Xk G S 1 - 5oi/t frame changes — /or in section l2J\ and 

9-dependent in section lKM — involve only phase changes on eigenstates. Therefore, for 
all t and for all \k) , 

\\\m)m)\-\k){k\\\ = \\m))m)\-\k)(k\\\ 

= \\m)m\-\k){k\\\. 

3. Robust ensemble transfer from \k) to \N — k — 1} 

In this section we consider a control protocol - chirped pulse - used by physicists to 
drive a system from the lowest eigenspace, spanned by |0), to the highest eigenspace, 
spanned by \N — 1), of the free Hamiltonian H given in fl2]). In fact we prove that 
a general (robust) class of control inputs transfers population from eigenstate \k) to 
eigenstate \N — k — 1), for all k, on an ensemble of systems with different values of 
parameters /x , • • • , (dipole moments) and A , . . . , Ajv-i (anharmonicities). 

The key requirements on the control are (i) to use a sufficiently chirped pulse - 
condition (b) in Theorem [1] — and (ii) to avoid all eigenvalue crossings — condition (c) 
in Theorem [H 

3.1. Transfer Theorem 

For k = 0, ... ,N - 1 let Af (s) = (k\H R (u(s))\k) = k(A k - cu(s)), the eigenvalues of 
H R (u(s)). 

Theorem 1. For given A > 0, fi max > Umin > 0, consider S an ensemble of systems 
of type fW\) with jij G [fi m in, Umax] for all j G Nq' 2 and Aj G [—A, A] for all j G Nq' 1 . 
Take controls A and uj with: 

(a) A and u G C 2 ([0, 1],R) 

(b) uj(0) and w(l) are such that, for all systems in S, 

A£(0) < ... < A^_!(0) and (13) 
A R (1)>...>A^ 1 (1) 

(c) A(0) = A(l) = and A(s) ^ for s G]0, 1[ 
Then 3 a constant C > such that, for all s > 0, 

sup || U £ (l)\k)(k\U £ (iy - \N - k - 1)(N - k - 1\ \\ < Ce . 
s 

The proof of this theorem is given in section EJ we there actually replace the simple 
condition (c) on A by a more general one: ^4(0) = A(l) = and A(s) ^ for all 
s G X w (5), where 

T u = {s G [0, 1] : Hr(uj(s)) is degenerate for some system G S} . (14) 
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The argument is based on the facts that the system approximately follows eigenstates 
of H(s) for small enough e (adiabatic theorem), eigenvalues of Hr are inverted between 
s = and s = 1 thanks to w(s) (chirping), and nonzero A(s) avoids all crossings for 
eigenvalues of H(s) such that e.g. the initial highest-energy level \N — 1) connects to 
the final highest-energy level |0) (see Lemma [1] in section [6]). Theorem [1] implies that for 
a given control satisfying the assumptions, taking e small enough allows to invert the 
state populations of a whole ensemble of systems featuring different parameter values. 
The control inputs only need to satisfy a few weak conditions and are therefore robust 
to many perturbations. These insensitivity properties of the adiabatic passage protocol 
have long been recognized by experimentalists. They commonly use the following type 
of control, see e.g. [T7] . 

Example 1. A function u satisfying the inequalities (T73p is e.g. u(s) = a(s — |), for 
a large enough positive a; such u is said to perform a frequency sweep. Except for the 
finite extension of time domain, such inputs are obtained by a "chirped" Gaussian laser 
pulse, which takes the form E(t) = Eq e~^ T2 e"^ e - ^' d( where k^O characterizes 
chirping. 

Theorem [1] still holds if inequality f JT3|) is replaced by 

A*(0) > . . . > A$_i(0) and A R (1) < . . . < A^(l) , 

i.e. the direction of the frequency sweep in Example []] can be inverted (taking a large 
enough negative a). However, for a given system, choosing one inequality over the other 
may allow to get a lower value for the constant hidden in the "order of magnitudes" 
result. This brings a mathematical foundation to the experimental observations made 
e.g. in [19]. 

3.2. Simulations 

We simulate system (TlOl) with a control satisfying assumptions (a), (b) and (c) of 
Theorem [U We consider a 4-level quantum ladder (so N = 4). We take e = 10 -2 , 
A ,...,A 3 G [-0.4,0.4] and ^0,^1,^2 e [/W, AW] = [1,5]. The control is u(s) = 
8(s — |) and A(s) = s(l — s), represented on Fig{T]a. FigfjTJb shows how the eigenvalues 
of H(s) (thick lines) avoid crossing For the illustrated random choice of detunings, the 
eigenvalues of Hr(u(s)) (thin lines) are very close to concurrent between s = 0.5 and 
s = 0.6. This poses no problem for the adiabatic transfer from \k) to \N — k — 1). 
The successful transfer is illustrated on Fig{T]d, which shows the squared norm of the 
projection of U e {l)\k) onto \p), for all pairs (\k), \p)) of eigenvectors of H ; this is 
equivalent to the squared norm of element on row p, column k of matrix U £ (l) that acts 
by left-multiplication on initial column- vectors, for U £ (l) expressed in basis (|0), . . . , |3)). 
Fig{T]c shows ensemble control on ten systems with different random values of A , . . . , A 3 
and no, 
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Figure 1. Control scheme transferring \k) to \N — k — 1). (a) control inputs A(s), 
uj{s). (b) s-dependent eigenvalues of H(s) (thick lines) and of Hr(iv(s)) (thin lines), 
(c) population on level |3) for 10 systems whose parameters ^0,^1^2 and Ai, A 2 , A 3 
were randomly picked respectively in [1, 5] and [—0.4, 0.4], and all starting at the initial 
state |0). (d) squared norm of the matrix elements of U £ (l), represented in shading 
from white (value 0) to black (value 1). 



4. Robust ensemble transfer from |/} to \p) 

In this section we propose a new robust control protocol to drive a system from the 
eigenspace (of free Hamiltonian H ) spanned by \l) to the eigenspace spanned by \p), 
for any given / and p in N^" 1 . The population transfer works on an ensemble of systems 
with different values of /zo, • • • , Hn-2 (dipole moments), and for a general class of inputs 
where zero- crossings of A(s) must be correlated with degeneracies of Hr(ou(s)); the latter 
depend on u(s) and (some of the) anharmonicities A , . . . , Ajv_i, which must hence be 
fixed. 

4-1. From |0) to any \p) 

For the sake of clarity, we start by giving sufficient conditions on A and u for the 
particular population transfer from |0) to arbitrary level \p). Section l4~2l generalizes the 
result to arbitrary initial state Consider the following assumptions: 

(Al) S is an ensemble of systems of type (TlTJT) with fij G [/i m i n , Umax] for all j G Nq' 2 , for 
some given fi ma x > Umin > 0, and with given sequence of detunings (A , . . . , Ajy-i), 
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such that the set {k(A k — v): k G 1 } contains at least A — 1 distinct values for 
any v G R; 

(A2) u is analytic and -£u)(s) > 7 > for all s G [0, 1]; 
(A3) cu(0) and cu(l) are such that ffTB"]) holds. 

For any m and n in N^" 1 with m < n, we note s(m, n) the unique timj§| where 
A^(s(m,n)) = \%(s(m,n)). 

As all systems in 5 have the same sequence of detunings, they feature the same 
eigenvalues A^, . . . , X^-i °f -^R an d hence the same set of s(m,n). The set of all 
s(m,n) equals X w defined in (JHJ), with dependence on particular system G 5 becoming 
irrelevant. The end of (Al) further implies that Hr has at most one pair of equal 
eigenvalues for any s G [0,1] i.e. (m,n) 7^ (j,k) implies s(m,n) ^ s(j,k), hence 
X w contains N(N - l)/2 distinct values. Further define 1$ = {si, . . . , s^v-i} C X w 
the N — 1 points where \q(s) = A^(s) for some n G N^ -1 , numbered such that 
si < s 2 < . . . < sjv-i- Thus, for each G Xq there exists a unique n G N^" 1 such that 
s k = s(0,n). 

The key requirements on the control to achieve population transfer from |0) to \p) 
are (i) to use a sufficiently chirped pulse frequency — condition (A3) — and (ii) to 
shape pulse amplitude in order to appropriately provoque — (c) in Theorem [2] — or 
avoid — (b),(d) in Theorem [2] — crossing of eigenvalues of H. 

Theorem 2. Consider S an ensemble of systems satisfying (Al) with a control u 
satisfying (A2) and (A3). Take p G {0,...,N — 1} and consider a control A with 
the following properties: 

(a) A is analytic over [0, 1] and A(0) = A(l) = 0. 

(b) A(s) ^ for all s G X^\X^. 

(c) A(s k ) = for all s k G 1q with k < N -p - 1. 

(d) A(s k ) ^ for all s k G X^ with k > N - p. 

Then 3 a constant C > such that, for all s > 0, 

sup || U £ (1)\0) (0\U £ (iy -\p}(p\ || < Cy/i. 
s 

The proof, given in section |6l shows that at eigenvalue crossing points the system 
adiabatically follows the eigenvector corresponding to the crossing branch. 

4-2. From any \l) to any \p) 

Under assumptions (Al) to (A3), we denote X£\(s) = {s(m,n) G T w : m = k, n > 
k and s(m,n) > s} and X^_(s) = {s(m, n) G X w : m < k, n = k and s(m,n) > s}, 
for any k G Nq' 1 . Further let q k ±(s) = inf(X^ ± (s)) and define g k ±(s) by s(k, g k+ (s)) = 

§ If assumptions (Al) to (A3) hold, then the existence and unicity of s(m,n) is ensured for all m and 
n > m: see Figj2]b or Fig(3jc 
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qk+(s) and s(<7fc_(s), k) = qk-(s) respectively. For p < N — I — 1, construct Xf p with the 
following algorithm. 

1: d := 0; x := 0; k := I; Xf p := 0; 

2: while d < N -I -p-l do 

3: while [Xfc_(x) ^ and < q k +(x) ] do 

4: fc := gk-(x)] x := g fc _(x); 

5: end while 

6: XI := Xf p U {q k+ (x)}; d := d + 1; x := g fc+ (x); 
7: end while 

The algorithm is verified to always successfully complete[jj]. For p > N — I — 1, we can 
define Xf p with a similar algorithm but where ' < ' is changed to ' > ' on line 2 and indices 
k--, k+ are switched. Then Xf p contains \N — I — p — 1| elements. 

Corollary 1. Consider S an ensemble of systems satisfying (Al) with a control oo 
satisfying (A2) and (A3). Take l,p in Nq' 1 and consider a control A with the following 
properties: 

(a) A is analytic over [0, 1] and A(0) = A(l) = 0. 

(b) A(s) =0 for allseX? p . 

(c) A(s) ^ for all s6r\I^. 

Then 3 a constant C > such that, for all e > 0, 

sup \\U%l)\l)(l\U%iy-\p)(p\\\ < CV~e. 
s 

Assumption (Al) ensures that each eigenvalue crossing / anti-crossing can be ad- 
dressed individually. This ensures that any transfer can be implemented in any situa- 
tion, but it is in general not necessary for a given system and transfer, as (simultaneous) 
crossings of some eigenvalue branches are irrelevant. The control proposed for Theorem 
[2] or Corollary [T] is just one amongst many possibilities of "eigenvalue crossing designs". 
Indeed, depending on (l,p) and on the particular arrangement of the s(m,n), one can 
find other subsets Ji v C X w such that taking A(s) = if and only if s e J~i p , permutes 
the eigenvalues in such a way that A;(l) = A^(l). The controls that we propose are 
optimal in the sense that they require a minimal number of pairwise crossings, that is of 
annihilations of A at accurate points. Variant annihilation subsets Ji v may be useful (i) 
to avoid some crossing points s(m, n) or eigenvalue branches (e.g. because corresponding 
A m or A n is poorly known, or because s(m,n) is close to some other point in T w ), (ii) 
to optimize adiabatic convergence as a function of e, or (iii) to simultaneously perform 
population transfers between several eigenstates, as we do in section [5j 



|| Indeed by construction, the cardinality of I^Ax) equals N — I — d — 1 (except during the update on 
line 6) and the cardinality of IP (x) decreases by one each time line 4 is applied; thus it is impossible 
to keep applying line 4 infinitely, and line 6 is always well-defined (that is IP + (x) ^ 0) for d < N — I — 1. 
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Another approach for transferring |Z) to \p) is to use A(s) Gaussian, i.e. without 
any annihilations, but reduce uj(s) to a specific range. Indeed, under the above 
assumptions, it is possible to choose w m i n and Umax such that /(A; — v) = p(A p — v) 
for some v G [w m j n , w m ax] and H R {v) is non-degenerate for all v G [w m j n , w m ax] \ {v}- 
Then taking uj(s) monotone between uj m ^ n and wmax just induces one avoided crossing 
that exchanges \l) and \p). Pictorially, this is like selecting a particular narrow vertical 
slice on FigJTJb. Depending on the specific system under study and whether it is 
experimentally easier to precisely modulate the amplitude or the phase of a field, one 
method may be more suitable than the other. A main advantage of our method is that, 
unlike the method proposed in [26], it can be extended to achieve any permutation of 
eigenstates as is shown in section |5j 

4-3. Simulations 




Figure 2. Control scheme transferring |0) to \2) control scheme; subplots analogue to 
FigfTJ except that A l5 A 2 , A 3 remain fixed for (c). A(s) vanishes at s = 0.25 so that 
Ao and Ai cross instead of avoiding crossing. 

As in section 13.21 we simulate (ITU]) for a 4-level quantum ladder (so N = 4) with 
p , fii, p 2 G [pmin, Umax] = [1,5]. We now take e = 1CT 3 and in accordance with the 
statement of Theorem |2] we fix the anharmonicities, taking Ai = —1, A2 = 0.3, A3 = 
(the value of Ao, multiplied by k = 0, is irrelevant). We target in particular a transfer 
from |0) to |2). The algorithm of section H~2l reduces to the simple case of Theorem |2l 
requesting a single zero of A(s) at s — inf {s(0, 1), s(0, 2), s(0, 3)} = s(0, 1) = 0.25 in 



Adiabatic passage and ensemble control of quantum systems 



12 



addition to A(0) = A(l) = 0. We take A(s) = s(l - s)(s - 0.25) and u(s) = 4(a - §), 
represented on Figj2ja. Figj2jb shows how the eigenvalues Afc(s) of i?(s) cross or not 
(thick lines); the eigenvalues of Hr(oj(s)) (thin lines) define points s(m, n) for our control 
design. Figj2jd confirms achievement of the intended result by showing the squared norm 
of components of matrix U £ (l) in basis (|0), . . . , |3)): we indeed have |(p|C/ e (l)|A;)| 2 ~ 1 
for (p,k) = (2,0) (other values incidental). FigfSJc illustrates ensemble control on ten 
systems with different random values of ^0)A*i;A*2- Since for this particular case the 
control only exploits precise crossing point s(0, 1) = 0.25, we might actually allow 
ensembles with different A 2 , A 3 . 

5. Robust ensemble permutation of populations 

In this section we describe the most general result of the paper, adiabatically transferring 
(|0)(0|,...,|iV-l)(iV-l|) to (\a(0))(a(0)\,...,\o-(N-l))(a(N- 1)|), where a is any 
permutation of N^" 1 . As in section HI the population permutation works on an ensemble 
of systems with different values of /^o, ■ • • , (dipole moments), and for a general 

class of inputs where zero-crossings of A(s) must be correlated with degeneracies of 
Hr(cu(s)); the latter depend on u(s) and require anharmonicities A , . . . , Ajv-i to be 
fixed and known. We prove existence of an appropriate control by recurrence on N. In 
fact this recurrence method can be used to design A(s), as we illustrate in section l5\2l 

5.1. Permutation theorem 

Theorem 3. Consider S an ensemble of systems satisfying (Al) with a control u 
satisfying (A2) and (A3). Take a any permutation o/Nq 1 - Then there exists a subset 
%a Q ^ for which, taking control A to satisfy 

(a) A analytic over [0, 1] and A(0) = A(l) = 0, 

(b) A(s) = for all s 6 T A , 

(c) A(s) ^ for all s e l u \ T A , 

implies: 3 a constant C > such that, for all s > 0, 

sup || f/ £ (l)|fc)(fc|f/ £ (l) t - \a(k))(a(k)\ \\ < C^fe. 
s 

Since the proof of this Theorem is constructive and necessary for the understanding 
of the example below, we present it here. 

Proof (of Theorem [3]) . The formal arguments (sup, adiabatic propagator) are presented 
in detail in the proof of Theorem Ul in section We focus on the construction of the 
control A(s) by following analytic eigenvalue branches of H{s). 
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The property is obvious for N = 2: either (cr(0), er(l)) = (1,0), which follows 
Theorem[J\just requiring A(s(0, 1)) ^ 0; or (cr(O), cr(l)) = (0, 1) , which follows Theorem 
d transferring |0) to \p) = |0) with one crossing^, i.e. just requiring A(s(0, 1)) = 0. 

Assume that we can achieve any permutation ofN^' 1 for N = K , and we are given 
a permutation a of Nq for N = K + 1 where a (I) = K and cr(K) = p. 

• If I = p = K , i.e. c{K) = K, then first build the remaining permutation on 
levels \0),...,\K — 1) by neglecting level \K). This uses the result for N = K; it 
just requires A(s) = for some s = s(m,n) and A(s) ^ for some other s = s(m,n), 
with m,n < K . Now take a particular such A(s) where in addition, A(s) = for all 
s G {s(m,K) : m G N^" 1 }. Then Xk(s), starting at Xk(0) = A^(0), exactly follows 
the same crossings as A^-(s) to end up as A^(l) = A^(l); the other levels remain 
unperturbed, so a is achieved. 

• If I 7^ K 7^ p, then first construct A(s) by applying the result of the preceding point to 
a, defined by 

a(l)=p; a(K) = K; a(k) = a{k) for all k {I, K} . 

A(s) performs the target permutation, except that K remains on K and I goes to p. From 
(E2P eigenvalue branch \k{s) necessarily crosses, at some s G {s(m,K) : m G N^" 1 }, 
the analytic eigenvalue branch that starts at A;(0) = Af (0) and ends at Xi(l) = A^(0). 
Define A(s) to have the same zeros as A(s) except that A(s) ^ 0. This just transforms 
the crossing ats into an anti- crossing, such that the analytic branch coming from A^(0) 
(resp. Xi(0)) now connects to the analytic branch going to A p (l) (resp. A^(l) y ). Thus 
A{s) achieves the target permutation a. □ 

Each "eigenvalue crossing design" choice X4 yields a particular permutation a% A - 
For N > 2, the number 2 N ^ N ~ 1 ^ 2 of possible X4 (i.e. subsets of X w ) is strictly larger 
than the number iV! of permutations. Thus there are still several X4 that yield the same 
a. Unlike in section HJ building A(s) as in the proof of Theorem [3] does not necessarily 
yield a minimal cardinality of X4 for given a. 

5.2. Example and simulations 

We first illustrate the control design by recurrence based on the proof of Theorem 
[3j Consider target permutation <r(0, 1,2,3) = (2,0,3,1). First we reduce it down 
to an elementary permutation. Start with K = N — 1 = 3 and note (l,p) = (2,1) 
because er(2) = K and cr(K) = 1; we thus define a(0, 1,/ = 2,3) = (2,0,p = 1,3) 
and impose A(s) = for s G {s(0, 3), s(l, 3), s(2, 3)} reducing the permutation to 
0,1,2. Then we take K = N — 1 = 2 and note (j,p) = (0,1) because cr(0) = K 
and a(K) = 1; we thus define a (J = 0, 1, 2, 3) = a(p = 1, 0, 2, 3) and impose A(s) = 
for s G {s(0, 2), s(l, 2)} reducing the permutation to 0,1. To implement a we need 
A(s(0, 1)) 7^ 0. Now we progressively move up to permutations on more levels, removing 



Indeed, {-Pa (i))-Pai(i)} = {|0)(0|, |1)(1|} then automatically implies transferring |1)(1| to P\ 1 (i) = 




Figure 3. Control scheme generating permutation er(0, 1,2,3) = (2,0,3,1) and 
simulation result. Subplots (a),(b),(d) analogue to FigO}. Subplot (c) shows the 
eigenvalues of Hr(uj(s)), used to design Xa (see text). The points s G 1a where 
A(s) = are marked on (a). 



one — at a time from our objects. The reader is encouraged to follow crossings/anti- 
crossings under the different controls by referring to FigJSJc, corresponding to our 
chirping choice oj(s) = 4(s — |). Under A the analytic branch from \l) = |0) to 
\p) = |1) meets the branch staying on \K) = |2) at s — s(l,2). We therefore impose 
A(s(l, 2)) 7^ unlike for A, and for the rest copy the requirements of A: A(s(0, 1)) ^ 0, 
A(s(0,2)) = 0. Now under A the branch from \l) = |2) to \p) = |1) crosses the 
branch staying on \K) = |3) at s — s(2, 3). We therefore get requirements for our 
actual control A by imposing A(s(2, 3)) 7^ unlike for A, for the rest copying the 
requirements of A, i.e. A(s) = for s G {s(0, 3), 3), s(0, 2)} and A(s) 7^ for 
s G {s(0, 1), s(l, 2)}. To satisfy these requirements, we take the polynomial control 
A(s) = s(l - s)(s - s(0,3))0 - s(l,3))(s - s(0,2)), represented on Figj3ja. FigfJJb 
shows how the eigenvalues of H(s) cross and anti-cross depending on whether A(s) 
vanishes or not. The squared norm components of U £ (l) resulting from a simulation 
of f fTUj) with this control and e = 10~ 3 are shown on Figj3)d on a white-to-black scale, 
confirming achievement of permutation cr(0, 1, 2, 3) = (2, 0, 3, 1). 

FigJUshows the same squared norm components of U £ (l) in gray-shades for 24 cases, 
corresponding to different control inputs A(s) designed for all 24 possible permutations 
of the set (0,1,2,3). The controls A(s) are built as the product of (i) a polynomial 
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Figure 4. Simulation results of (fTU| for 24 different controls A(s) following Thcorcm[3] 
to achieve each one of the 24 permutations of (0, 1,2,3) with adiabatic passage. Shading 
represents squared norm of elements of matrix U 6 (l) expressed in basis |0), . . . , |3), 
from white (value 0) to black (value 1). In other words, each subplot may be read as 
a 4 x 4 matrix where the black patches are ones and the white patches are zeros; gray 
patches indicate intermediate values, reflecting that the unitary propagator obtained 
by integrating (|10p is not exactly a permutation matrix for the finite e = 10~ 3 . 



vanishing on Ia U {0, 1} and only there, and (ii) a set of functions (1 + g(s — s(m, n))), 
with g(s — s(m,n)) Gaussians centered on all s 6 I u \ Z^; the role of the latter is to 
amplify A(s) in the vicinity of intended "anti-crossings", improving convergence of the 
adiabatic limit as a function of e. FigH] corresponds to the choice e = 10 -3 . 



6. Proofs 



In this section we give the proofs of all the formal results presented in previous sections. 



6.1. Proof of Theorem^ 

We start the proof by recalling the following result |27j . 

Lemma 1. Let be a real tridiagonal and symmetric N x N matrix defined by 

N-l N-2 

D N = J2 a k\k)(k\ + c k(\k)(k + 1| + |* + l)(k\) (15) 

fc=0 fc=0 

in some orthonormal basis (|0), . . . , \N — 1)). If Ck ^ for all k G Nq' 2 , then Dn is 
non degenerate. 

Proof (of Lemma [T]). Denote Q n the characteristic polynomial of D n , which is defined 
as l[T5\) with N replaced by n, for alln G {1, . . . , iV}. The sequence of polynomials (Q n ) n 
verifies the following recurrence relation: for n > 2, 

Q n {x) = {x- a n _i)<5 n _i(x) - (C n _ 2 ) 2 Qn-2(x) , (16) 
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with Qo(x) = 1 and Q±(x) = x — a . According to Favard's Theorem a sequence 
verifying j[Tb)) where (c n _ 2 ) 2 > for all n, is a sequence of orthogonal polynomials. 
Furthermore, from J2E. Theorem 3.3.1], every polynomial Q n in a sequence of orthogonal 
polynomials has n real and distinct zeros; this is in particular true for n = N, therefore 
Dn is non degenerate. □ 

Proof (of Theorem [T]) . We prove the result for any single system in S and conclude 
that it remains true for the sup over S. Indeed, the application 

(fj, , . . . , fiN-2, A , • • • , A N _i) 

-> || U £ (l)\k)(k\U e (iy -\N-k- 1)(N — k — l\ || 

reaches its sup over the allowed compact space since the state of a (sufficiently regular) 
dynamical system at a finite time depends continuously on system parameters (see 
e.g. WB, theorem 3.5]). The proof for one system is in two steps: first we prove 
that the hypotheses of the adiabatic theorem with gap condition are verified, then we 
apply the theorem to compute the image at s = 1 of initial projector \k)(k\ in adiabatic 
approximation. 

Step 1: By hypothesis (a), we have H(s) G C 2 ([0, 1],"Hjv) an d therefore continuous 
over [0, 1]. From IUffl, section 11.5.2], it is then possible to find N continuous functions 
Ao(s), . . . , \n-i(s) such that Ao(s) < . . . < Xn-i(s) for all s G [0, 1] are the eigenvalues 
of H(s). In terms of associated eigenspace projections, note that {P\ k ( s ) ■ k G N^" 1 } = 
{|A;)(A;| : k G N^" 1 } every time A(s) = 0, by unicity of the spectral decomposition of a 
non degenerate matrix. However, the pairwise correspondence depends on the value of 
id(s). In particular, by hypotheses (b) and (c), 

P Xk[0) = \k)(k\ and P Xk(1) = \N-k-l)(N-k-l\ (17) 

for all k. For a given s G [0, 1], 

• either A(s) ^ 0, then H(s) has N distinct eigenvalues according to LemmaUk 

• or A(s) = 0, then H(s) = H R (u(s)) and it must have N distinct eigenvalues by 
hypothesis (c). Hence, 

A (s) < < Ajv-i(s) for all s G [0, 1] . (18) 

Then by continuity over the compact [0,1], there exists 5 > such that Afe(s) + 5 < 
Afc+i(s) for all k G N^" 2 and for all s G [0, 1]: each Afc(s) is at all times surrounded by 
a "spectral gap" of amplitude 5 in which there is no other eigenvalue. We can therefore 
apply the adiabatic theorem with gap condition (see \12\ Theorem 2.2]) to eigenvalue 
Afc(s), for any particular k G N^" 1 , as is done in the following. 

Step 2: The adiabatic theorem ensures that P Xk{s) G C 2 ([0, 1],H N )- Define the 

"adiabatic Hamiltonian" 

H a>k (s) = H(s) - isP Xk{s) ±-P Xk{s) - iePt(s)-^Pt(s) (19) 
where Px k r s ) = I ~ ^A fc (s); an d the "adiabatic propagator" U e ak which verifies, for all 

se [0,1] 

is ^ Ul k (s) = H a!k (s) Ul k {s) with UIM = I . (20) 
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One verifies that this construction ensures 

UlM P MO)U £ aM = Px k {s) (21) 

for all s G [0, 1]. The adiabatic theorem states the existence of a constant C\ > such 
that 

||C/ £ ( S )-^, fc (s)||<C l£ for all se[0,l], 
in particular for s — 1. This implies 

||^ £ (l)|A;)(A;|^(l)t-^ fc (l)|A;){A;|C/^(l)t|| 

< ll(^(i) -^(i))!*)^!^!! 

+ \\u £ a Ai)\k)(k\(u £ (i)-ui k (i)y\\ 

<\\U £ (1)-U^ k (l)\\ \\\k)(k\\\ (11^(1)11 + 1^(1)11) 

< de • 1 ■ 2VN 

since \\U\\ = Vtr WU = vtrJ /or any unitary matrix U. Combining this with (Tty, $2~1\) 



yields the result, where C = 2C\\fN . □ 
6\£. Proof of Theorem^ and corollary^ 

We start by proving a Lemma about the behavior of time- dependent eigenvalues crossing 
each other. 

Lemma 2. Assume that H[s) as defined in ( TO]) depends analytically on the real 
parameter s on an interval Icl, with 4-w(s) > 7 > for all s E T. Suppose that 
Hr(uj(s)) is non degenerate on X except for a simple degeneracy atsET, i.e. Hr(oj(s)) 
has N — 1 distinct eigenvalues and Hr(uj(s)) has N distinct eigenvalues for s G T\ {s}. 
If A(s) = 0, then: 

(a) There exist N unique functions Xq, . . . , Aat_i analytic over T, with Xq(s) < . . . < 
Ajv-i(s) for all s < s, and such that {Xq(s), . . . , Ajv-i(s)} the eigenvalues of 
H(s) for all s G X. 

(b) Let k be such that Xk{s) = Xk+i(s). Then for all s > s we have 

X (s) < ... < X k+1 (s) < A fe (s) < . . . < Ajv-i(s) . 

Proof (of Lemma [2]) . Point (a) is a direct consequence of 12^ Theorem 6.1]. The order 
of the analytic eigenvalues is obviously preserved over time intervals where H(s) is non 
degenerate; by LemmaUl these intervals are {s < s} and {s > s}. The issue is what 
happens at s = s. In the following, we show that X' k (s) 7^ X' k+1 (s). Since the eigenvalues 
are analytic and X k (s) = X k+ i(s), a Taylor expansion then yields the conclusion of (b). 

We lead calculations similar to those of fTB. section XVI. II. 8]. According to 
section II. 6. 2], since H is analytic over T and H(s) G Hn for all s G X, there exist 
rank one orthogonal spectral projections P\ ( s ), ■ ■ ■ , -Pa jv _ 1 (s) which are analytic over X. 
Computing the derivative of 

H(s)P Xkis) = X k (s)P Xk(s) (22) 
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with respect to s at s = s, we get 

H'(s)P Xk{s) + H(s)P' Xk{E) = X' k (s)P Xkis) + h(s)P' Xk{ - s) ■ 

Multiplying the last equation by {P\ k (s) + P\ k+1 (s)) from the left, using $2E) and the 
fact that P\ k and P\ k+1 are two orthogonal projectors (P Xk = P\ k , P\ k 1 = P\ k+1 
and P\ k Px k+1 = 0), we get (P Afc(g) + P Xk+l{s) )H'(s)P Xk (s) = K(s) P \ k (s)- Noting that 
P\ k (s) = (P\ k (s) + Px k+1 (s))P\ k (s), we get 

(PXk(s) + P\ k+1 (s))H'(s)(P Xk{ s) + P\ k+1 (s))P\ k (s) 

= K(*)Px k & ■ 

The analog holds with k and k + 1 switched. This implies that {A^(s), A / fc+1 (s)} are the 
eigenvalues of the 2x2 matrix obtained by restricting operator H' '(s) to the column space 
of (P\ k (s) + P\ k+1 (s))- Since A(s) = we have H(s) = H R (u(s)). Denoting \m) and 
\n) the two eigenvectors of Hr corresponding to eigenvalue Xk{s) = \k+i{s), we have 
P\ k {s) + P\ k+1 (s) = \m)(m\ + \n)(n\. Defining 



(H'(s)) mn 



(m\H'(s)\m) (m\H'(s)\n) 
(n\H'(s)\m) (n\H'(s)\n) 



and computing 



H f (s) = u'(s)^-H R (v) \v=u{s) +A'(s)Hi , we get 

{H {S)) - ~ { A>(s)» mn -W(S) J (23) 

where [i mn = (m\Hi\n) . Thus \i mn = if \m — n\ > 1 and fi mn ^ if \m — n\ — 1. 
In both cases, since uj'{s) ^ and m ^ n, the matrix in ^2^1 has 2 real and distinct 
eigenvalues, corresponding to A' fc (s) 7^ A' fc+1 (s). □ 

Proof (of Theorem |2]) . Taking A(s) = at some points where Hr is degenerate means 
that eigenvalues of H(s) will not remain distinct at those points. We therefore use 
the adiabatic theorem without spectral gap condition, see 11^ corollary 2.5]. Like for 
Theorem^ we prove the result for one system £ S and conclude the result for the sup. 
The proof is again in two steps. First we state how the adiabatic theorem can be applied; 
then we compute the image at s = 1 of initial state \k)(k\ in adiabatic approximation. 

Step 1: Since H is Hermitian, analytic over [0, 1] and simply degenerate at isolated 
points, we can apply Lemma W(a) repeatedly to conclude that there is a unique set 
of functions Ao, . . . , Ajv-i analytic over T, with Ao(0) < ... < Ayv_i(0), and such 
that {Ao(s), . . . , Ajv-i(s)} are the eigenvalues of H(s) for all s £ X. Moreover, 
according to 12$, section 11.6.2], there is a unique set of associated rank-one projectors 
P\ (s), ■ ■ ■ 1 -Pajv-i(s) which are analytic over T. In particular, given assumption (A3) 
and as H(s) = H R (u(s)) for s £ {0, 1}, we have (A fc (0), P Afc(0) ) = (Af (0), for 
all k and {(A (l), Pa o( i)), • • • , (Aiv-i(l), = {(A R (1), |0)<0|), . . . , (A^l), |JV- 
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l)(iV — 1|)}. Note however that, unlike for Theorem [7J the pairwise correspondence 
between elements of the latter sets is not obvious a priori, because here eigenvalues of 
H(s) do not remain distinct on [0, 1]. A second difficulty is to assess how the system's 
state evolves when eigenvalues become degenerate. This second part is answered by 
the adiabatic theorem witout gap condition. Introduce, as in TheoremUl the adiabatic 
Hamiltonian H a Q and adiabatic propagator U^ , given by / TTPj) and respectively with 
k = 0. Then by construction ^, (l)|0><0|f/ a yi)t = U^l) P Xo{0) Ul (iy = ^A (i)- The 
adiabatic theorem states that 3C such that 

\\U%s)\k)(k\U^ -E£ ( S )|*)<*|E£o( a )t || < C^e (24) 

for all \k) G {|0}, . . . , \N — 1}}. Thus the actual system adiabatically follows the analytic 
P\ k ( s ), from P Ao( o) = |0)(0| up to P A() (i) in particular. 

Step 2: We now compute -Px (i)- Define a small interval I mn = [t^^t^J C [0,1] 
around each point s(m,n) such that all I mn are disjoint. If A(s(m,n)) ^ 0, then as 
shown in Theorem [JJ, H{s) is non degenerate for all s G I mn , such that for any j, k 
with Xj(r^ n ) < X k (r^ n ) we have Aj(r^ n ) < A fc (r£ n ). On the other hand, if we take 
A(s(m,n)) = 0, then two eigenvalues intersect at s = s(m,n) and the analytic branches 
cross so that their order changes as stated in Lemma\^(b) . To avoid separate treatment 
of limit cases, we define sq = and sjy = 1- Now by construction: 

• Aj(0) = Xf (0) for all j G Nq' 1 . 

• For k G Nf, Xq(s) is the kft 1 smallest eigenvalue of H R (u(s)) when s G (sfc_i,Sfc). 

• As long as A(so) = . . . = A(sk-i) = 0, that is for k < N — p, Xo(s) follows the 
same crossings as X{j(s); therefore it is the k smallest eigenvalue of H(s) when 
s G (sfc_i, s fe ). 

• For s > sjv-p-i, we have A(s) ^ so the Xk(s) keep the same order, i.e. Xq(s) 
remains the (N — p) smallest eigenvalue of H(s). 

• In particular for s = 1, from ([73)) we identify Ao(l) = X^ } _^ N _ p - ) (l) = X p (l), such 
that -Pxo(i) = \p)(p\ by uniqueness of the spectral decomposition. □ 

Remark 2. To apply the adiabatic Theorem \T% corollary 2.5], it is sufficient to have 
P\ ( s ) £ C 2 ([0, 1],T-Ln)- However, a condition like H(s) G C 2 ([0, does not ensure 

the existence of P\ ( s ), ■ ■ ■ , -P\jv-iW e ^ 2 ([0' M^n), see 12^, example 5.3]. It is only for 
analytic H(s) that we can guarantee analytic P\ ( s ), which then in particular belongs to 
C 2 ([0,1],^). 

Proof (of corollary [T]) . The arguments are the same as in the proof of Theorem [H We 
concentrate on tracking the analytic eigenvalue branches Ao(s), . . . , Ajv-i(s) of H{s) to 
establish their pairwise correspondence with eigenvalues A^(s), . . . , A^-_ 1 (s) of H r (oj(s)) 
at s = 1. We prove the result for p < N — I — 1; the case p > N — I — 1 is treated 
similarly, while p = N — I — 1, implying Tf p = 0, is the case covered by TheoremUl 
Denote s± < . . . < Sjv_;_ p _i the elements ofTf p , and s$ = 0, sat_;_ p = 1. 

The algorithm constructs Tf p such that the {I + d) and (I + d + 1) smallest 
eigenvalues of H R (u(s)) become equal at Sd, for each d G ^~ l ~ p ~ l . Taking A(sd) = 



Adiabatic passage and ensemble control of quantum systems 



20 



implies H(sa) = Hn(oj(sd)) so the same eigenvalue equalities hold for H(s) at s = Sd- 
Moreover from point (c) and Lemma [J\ all eigenvalues of H(s) remain distinct for 
s (jL Xf . Therefore the analytic eigenvalue branch Aj(s), starting with A;(0) = Af (0), 
exactly evolves through crossings at s\, . . . , s^-i-p-i such that it is the (I + d + 1)^ 
smallest eigenvalue of H(s) for s E (sd, s«m-i)- In particular, A/(l) is the (N — pj^ 1 
smallest eigenvalue of H(l) = Hr(uj(1)), which from (A2) means A/(l) = A^(l) such 
thatP Xl{1) = \p)(p\. □ 



7. Summary and discussion 

This paper shows how adiabatic passage can be applied to a quantum ladder system 
to achieve permutations of populations on the ladder levels with a single laser pulse. 
We explicitly propose control inputs whose precise functional dependence on time is not 
important as long as they satisfy a few key features, most notably annihilation or not at 
specific times. This makes our strategy robust against multiplicative input disturbances. 
Another important advantage of our adiabatic strategy is its ability to simultaneously 
control an ensemble of systems with different dipole moment values. 

Theorems in the present paper provide a proof of concept in idealized situations. 
Several practical issues deserve a more quantitative investigation in future work. 
Probably the most important aspect is to characterize precision of the adiabatic 
approximation as a function of e. Indeed, for small e the actual control time t — | 
gets long; this further implies that, at constant power A 2 (s), the energy given to the 
system gets large. Beyond performance requirements, this also invalidates our model at 
infinitesimal e (e.g. regarding finite lifetime of the levels). Although orders of magnitude 
are given for the adiabatic limit, variations in the proportionality constant can lead to 
significant discrepancies. Investigating them, as well as "optimal paths" minimizing 
non-adiabatic losses [30], could yield guidelines for choosing amongst several possible 
"eigenvalue crossing designs" . Both precision of adiabatic approximation and modeling 
assumptions (e.g. RWA) also limit the range of "ensemble" properties in practice. 

It may appear surprising at first sight that two different evolutions are selected 
just by taking A(s) = or A(s) ^ at a precise instant s. The elucidation is that 
this dichotomy only holds at the limit e — > + . For a given e, the larger \A\ in the 
neighborhood of s = s, the more the evolution differentiates from the A(s) = case. 
Nevertheless, for small e, the relevant neighborhood around s for selecting population 
transfer or not indeed gets small (from there experimentalists' denomination "rapid 
adiabatic passage"). Our scheme might therefore allow selective population permutation 
as a function of {Ao, . . . , Ajy-i} on an ensemble of systems, in a scheme resembling 
resonance selection. Take A(s) = for s G I u of a nominal system. If a system has 
detunings very close to nominal, then two of its Af (s) cross at a point s close to s, 
where A(s) ~ 0, such that for moderate e its final state will be close to the adiabatic 
result of the nominal system with A(s) = 0. If a system has detunings more different 
from nominal, then all its crossings of A^(s) occur at points where A significantly differs 
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from zero, and with moderate e its final state will be closer to the adiabatic result of the 
nominal system with A(s) ^ 0. A quantitative statement of this idea requires further 
investigation. 
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